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Abstract 

We study the Mott transition in a two-dimensional lattice spinless fermion 
model with nearest neighbors density-density interactions. By means of a 
two-dimensional Jordan- Wigner transformation, the model is mapped onto 
the lattice XXZ spin model, which is shown to possess a Quantum Group 
symmetry as a consequence of a recently found solution of the Zamolodchikov 
Tetrahedron Equation. A projection (from three to two space-time dimen- 
sions) property of the solution is used to identify the symmetry of the model 
at the Mott critical point as U q (sl{2)) (g> U q (sl(2)), with deformation param- 
eter q = — 1. Based on this result, the low-energy Effective Field theory for 
the model is obtained and shown to be a lattice double Chern-Simons theory 
with coupling constant k = 1 (with the standard normalization). By further 
employing the Effective Filed Theory methods, we show that the Mott transi- 
tion that arises is of topological nature, with vortices in an antiferromagnetic 
array and matter currents characterized by a d-density wave order parameter. 
We also analyze the behavior of the system upon weak coupling, and conclude 
that it undergoes a quantum gas-liquid transition which belongs to the Ising 
universality class. 
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1 Introduction 



In spite of many substantial advances in recent years, the study of some problems 
in the physics of strongly correlated electrons continues to provide stimulating chal- 
lenges. Among the several reasons for this, we could mention the still incomplete 
understanding of all of the properties of strongly correlated electron systems and the 
lack of many reliable techniques to study them. One of these paradigmatic prob- 
lems is the Mott transition, loosely defined as a metal-insulator transition driven 
by correlations. As early as in 1939, Mott argued that if the electronic density in 
metallic systems was lowered enough, the Coulomb repulsion would dominate over 
the kinetic energy and the system could undergo a transition to an insulating state 

[!]■ 

On the one hand, many widely employed theoretical descriptions of the Mott 
transition are based on the study of the microscopic dynamics of the electronic 
system: one starts by writing down a model Hamiltonian for the electrons, sometimes 
coupled to external fields o to some other degrees of freedom, and then tries to solve 
it within some approximation or through the use of numerical methods. However, 
numerical methods like exact diagonalization are restricted to small clusters, and if 
the interaction among electrons is strong enough, the typical approximation schemes 
based on the resummation of some class of Feynman diagrams are not completely 
reliable. For these reasons, it is a difficult task to find solutions displaying the 
Mott transition, even in the simplest models, like the Hubbard one. As of today, 
there exist two successful approaches for study this transition based on microscopic 
dynamics: one is the dynamical mean field theory (DMFT) method [2], valid in 
the limit of infinite spatial dimensions, which neglects the spatial correlations; the 
second approach consists in finding analytic expressions for physical observables 
in integrable models (mainly in one-dimensional systems) by means of the Bethe 
Ansatz or bosonization methods. 

On the other hand, systems like conventional superconductors or quantum Hall 
systems have universal properties that are well described by field theories which do 
not deal with the microscopic degrees of freedom, but rather with fields representing 
effective degrees of freedom. These two cases are classic examples of the more general 
framework of the Effective Field Theories (EFT) approach, which has its roots in 
Landau's ideas for condensed matter systems and which is widely and successfully 
used in high energy physics. This approach can be considered as 'way of thinking' 
which emphasizes the symmetries of the systems and that naturally incorporates 
Wilson's renormalization group ideas. In this approach, the study of a system starts 
by wisely choosing the effective degree of freedom, which are the relevant ones at a 
given energy scale, and then one proceeds to write down the most general second 
quantized action compatible with the characteristic symmetries of these degrees 
of freedom, retaining only the marginal and relevant terms, i.e., terms that are 
non-decreasing in the low energy (long-distance) limit [3]. During the last years, 
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several exotic states that contain droplets of approximately constant density have 
been found experimentally in electronic systems considered to be close to the Mott 
transition. These states that may seem to be surprising and difficult to explain from 
the point of view of the free electrons, are good candidates to be understood following 
the EFT approach. For example, the appearance of effective gauge forces arising 
from the dynamics, and which have not been included in the microscopic electron 
Hamiltonian can be properly taken into account within the EFT framework[4]. 

The goal of this article is to construct and consider an EFT for a two-dimensional 
square lattice system which displays the Mott transition, which implies that we 
shall focus our attention on the symmetry aspects of this transition. Specifically, 
we will consider a simple model of electrons with nearest neighbors density-density 
interaction which has also been previously studied, with the goal of identifying its 
effective degrees of freedom and their characteristic symmetries. Since the scope 
of this article is to make it readable to both Condensed Matter and Field Theory 
physicists, we shall also review (without pretending to be exhaustive) some basic 
aspects of bosonization, Conformal Field Theory, and Integrable Models. 

The paper is organized as follows: in Section 2, we present the model of strongly 
coupled fermions on a (two-dimensional) square lattice that we shall consider in the 
paper. We review some known properties of the corresponding one-dimensional ver- 
sion of this model and we also discuss the relationship among different approaches 
for treating the one-dimensional Mott transition. We also apply a two-dimensional 
bosonization prescription as considered in [5] for the two-dimensional fermion model. 
After reviewing some basic properties of integrability in statistical mechanics models, 
we discuss in the Section 3 the integrability of the specific two-dimensional fermion 
model considered in this paper. We show that in the strong coupling regime, the 
system defined by the ground and low-lying states of the model satisfies the Zamolod- 
chikov Tetrahedron Equation, and is characterized by a novel family of solution to 
the Tetrahedron Equation recently found by Bazhanov et. al. [6] [7]. We review 
these solutions for the sake of completeness and discuss the three-dimensional struc- 
ture of an underlying Quantum Group algebraic structure. This analysis allows us 
to identify the symmetry of the model at the Mott transition point as given by the 
Quantum Group U q (sl(2)) <8> U q (sl(2)). The identification of the symmetry and the 
corresponding effective degrees of freedom allows us to write down the EFT for the 
model, which is done in the Section 4. With the EFT at hand, we then analyze the 
order parameter and the universality class of the transition. We find that it is given 
by a Kosterlitz- Thouless type transition, with vortices in an anti-ferromagnetic ar- 
ray. We also discuss how this order is modified by doping, and that this procedure 
induces an Ising-like phase transition. Finally, we present our Conclusions. 
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2 The two-dimensional fermion lattice model 



In this Section we introduce the model that we shall be considering throughout the 
paper. Let us consider a spinless fermions system with nearest neighbors interaction 
on a square lattice, with Hamiltonian 

H 2d = -- Yy ^{x + ae^)e iA ^{x) + h.c. ] + U ^ p{x)p{x + aej , (2.1) 

X,fl X,fJ, 

where ip(x) is the fermionic field, x labels the lattice sites and are the unit 
lattice vectors pointing to the nearest neighbors of a given site, a is the lattice 
spacing , t is the hopping parameter, U is the (constant) Coulomb potential, p(x) 
is the charge density (normal-ordered with respect to the half-filling ground state), 
p(x) = [ : ^(x)i>(x) : -1/2] and is an Abelian statistical gauge field defined on 
the links of the lattice. 



2.1 Review of the one-dimensional model 

In order to proceed in our study of this model, we first would like to review the 
physics of the one-dimensional model analog of (2.1), given by the Hamiltonian 

Hid = ~\ Y}^( x + a )^( x ) + h - c -] + U J2p(x)p(x + a) , (2.2) 

~ X X 

were the sums are taken over the lattice sites. Note that the gauge field is unimpor- 
tant in this case, as it should, given that there are no statistical Gauge fields in one 
spatial dimension (see, e.g., [8]). This model has an interesting history that begins 
with the work of Luther and Peschel [9] and which has later on been studied in detail 
by several authors, including Shankar[10]. The Mott transition in one-dimensional 
systems has been discussed not only in the context of Hamiltonian models like (2.2), 
but also within the scope of Luttinger liquids. In the following, we review and relate 
both of these approaches with the scope of setting up a framework suitable for fur- 
ther generalizations and for finding the effective degrees of freedom for the simplest 
one-dimensional case. 

As it is well-known, the model (2.2) can be mapped onto the (one-dimensional) 
XX Z model through the Jordan- Wigner transformation: 

S + (i)=4exp(znJ2 1>fyi) , (2-3) 

i<j 

S-(t) = exp(-t7rJ2^ J )^, (2.4) 

i<j 

5,(0=^(0^(0-1/2, (2-5) 
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where % labels the lattice sites. The Hamiltonian goes onto 

Hit = £ [ -OS^i + + A5*(t)5*(< + 1) ] . (2.6) 

i 

In the thermodynamical limit, when the total number of sites is even, this Hamil- 
tonian corresponds to a the system at half-filling, which has the property: 

H XXZ (A) = -H xxz (-A) . (2.7) 

Moreover, the XX Z model is integrable and its spectrum and other analytic prop- 
erties have been found in [11] using the Bethe Ansatz. Its ground state energy 
is: 



\ cosh A - | sinh A [A + 2A £ n (l + e 2 ^)- 1 ] if A = cosh A > 1 

E Q = { 1/4 -In 2 ifA = l (2.8) 

j cos p — sin 2 p dx/[2 cosh7rx(cosh 2px — cos p)} if A = cos p < 1 

Note that Eg is an analytic function of A in the range 1 < A < oo, so that the 
singularity at A = 1 signals a phase transition. This fact has been used in [10] to 
show that the transition at at A = 1 is identified as a Mott one. The argument 
relies on the duality between two opposite regimes for the system, ranging from the 
insulator behavior for A — > oo to the metallic one for A = 0. Moreover, it is known 
that for A = 1 + e the spins are in a Neel state, and therefore the system must be 
in a charge density wave (CDW) state. 

An alternative description of the system (2.2) is given by bosonization of its 
fermionic degrees of freedom. For A < 1, the action of the system is given by 

S = -L I dzdz d z <P{z) d- z 4>(z) , (2.9) 
Aix J 

where g is a self-coupling parameter and we have defined complex space-time coor- 
dinates z = x + it and z = x — it, and the normal ordered charge density is given by 
p(z) = id z (f)(z) . The model has effective degrees of freedom which are bosonic fields 
representing charge density waves. Eq. (2.9) defines a conformal field theory (CFT) 
whose energy-momentum tensor has holomorphic and anti-holomorphic components 
given by: 

T(z) = -gd z (t>{z) (2.10) 
f(z) = -gd-J(z) . (2.11) 

Moreover, the Fourier modes of the fields, defined by 

T(z) = Y,L n z- n -\ (2.12) 

n 

p{z) = Y J P n z- n - 1 , (2.13) 
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satisfy the following chiral algebra: 



[ L n , L m ] = (n-m) L n+m + — 5 n+mfi (n 3 - n) (2.14) 

[ Pm Prn ] = H> $n+m,o (2-15) 

[ L m ,p n } = -m p n+m . (2.16) 

These three lines define a current algebra: the first line is the Virasoro algebra for 
the generators L n with central charge c — 1. The second one is the u(l) current 
(or Kac- Moody) algebra for the charge modes. The third is required for consistency 

among the other two. The u(l) current algebra is used to define the Luttinger 
model, which in Hamiltonian form is usually written as: 



/ dx L^(vrn(:r)) 2 + ^ (d^) 2 } , (2.17) 
!7T J L K J 



H = — 

2tt 

where u , K are called the Luttinger parameters. It is straightforward to regain the 
action (2.9) with coupling constant g = K starting from the Hamiltonian (2.17), by 
transforming the arguments of the fields into imaginary time. The parameters of 
the model in the different representations are related by (for a detailed discussion 
see [12]): 

K = - -— (2.18) 

2[vr - arccos(A)] v 1 

Hence, the Mott transition in the one-dimensional lattice fermion system (2.2) is 
characterized byA = lor.K' = g = l/2. The Luther-Emery transformation [13] 
allows us to rewrite the Hamiltonian (2.17) as H — H + H 1: where: 

H = v Jp &i(p)Mp) ~ 4>Hp)Mp)} dp (2.19) 
H x = ™ sinh(2£)[ [ (2 P+P - +fiJ2 : Pa(p)Pa(-p) -)dp } , (2.20) 

a=± 

where v = -u(cosh2# + fi sinh26 l ), exp(2#) = 1/{2K) and fi is an arbitrary constant. 
Note that for the Luther-Emery line, which coincides with Mott transition, i.e., 
K = 1/2, one has that Hi vanishes so that the theory consists of two free decoupled 
chiral fermions. Therefore, we identify the effective degrees of freedom of the theory 
at the Mott transition as two free fermionic currents of opposite chirality. 

Next we would like to discuss the characteristic symmetry of the degrees of free- 
dom and, therefore, of the system at Mott transition point. In order to do this, 
it would be more convenient to switch to the spin representation {i.e., the XX Z 
model). Note that the relation (2.7) allows us to write 



TTld . 

XXZ / j 



( QX QX I QU Q \ i,Q Q ) QZ QZ 



(2.21) 
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where q = expry, 7 = cos A i.e., A = — (q + q~ 1 )/2. For q = — 1, we have the 
isotropic anti-ferromagnetic spin chain, which is a critical system with an explicit 
SU(2) symmetry. 

It is well-known that the corresponding low energy effective field theory is the 
Wess-Zumino-Witten (WZW) model with coupling constant k = 1 [14] and action: 

Swzw = fd 2 x Trid.gd^g- 1 ] + / d^xTr^g'^gg-^gg^d^] , 
4tt J I2n J b{d) 

(2.22) 

where g(z) is a field in the group manifold of SU(2) (i.e., g is a S , f/(2)-valued 
matrix). The first integral in (2.22) is defined over a compactified two-dimensional 
domain D and the second is done over a three-dimensional ball with boundary D. 
For the sake of completeness, we will sketch here the procedure leading to this EFT. 
Following [14], one makes a transformation from the spin variables Si to a fermionic 
system ijjf that preserves the SU(2) symmetry defined by Sj = 1/2 ^ a ~ ■i/j^aip? 
, where o = (a x ,a y ,a z ) and <jj are the Pauli matrices. The standard commutation 
relations reproduce the correct spin commutators, however the Hilbert space of the 
fermion system is too large and one must restrict it by projecting out the states 
with one particle by site 

(^TOn)a = 1 • (2.23) 

In the low energy regime, the only excitations that should be taken into account are 
localized around the two Fermi points of the one-dimensional Fermi surface. One 
transforms to a set of new fermionic degrees of freedom (we will label the lattice 
sites by the integer n to avoid confusion): 

Mn) = \i n ^a L (n ± 1/2) + (-1)^(71 ± 1/2)] . (2.24) 

A final redefinition of variables takes us to current operators: 

J = 1 : (^iTMa : J l = i: (if>l) a A*l>L) a ■■ (2.25) 
G = + *PUl ^ L a^ R : + : ^ R a^ L : , (2.26) 

where R and L denote the left and right chiral components. The spin operators and 
the constraint (2.23) become: 

J + J = G = (2.27) 
S l /a = J l + J l + (-l) n G l . (2.28) 

In the continuum limit, the Hamiltonian (2.21) of the spin chain becomes: 

H = I J d 2 x [ f{x)J l {x) + J l {x)J i {x) + 2 fj l {x) ] , (2.29) 

where the x variable is the continuum limit of the lattice position. The last term in 
(2.29) is irrelevant in the renormalization group sense and the effective Hamiltonian 
becomes the Hamiltonian of the WZW model. 
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2.2 Bosonization of the two-dimensional fermion model 

We now turn our attention back to the two-dimensional fermionic model on the 
square lattice with Hamiltonian (2.1). It can be bosonized as discussed, e.g. in 
[15] and [16], by applying a two-dimensional Jordan- Wigner transformation. Let 
us first consider the case U — 0. As it is known from the one-dimensional case, 
the Jordan- Wigner transformation owns its existence to a natural ordering of the 
particles along the line. This ordering is lost in on two-dimensional lattice, but 
the mapping could still be defined by adding extra degrees of freedom, in the form 
of attached statistical fluxes to the particles, i.e., by introducing branch-cuts on 
the otherwise analytic fermionic field operators. Equivalently, one considers the 
Hamiltonian (2.1) with the additional Gauss law constrain: 

p(x) - 9B(r) = , (2.30) 

where p(x) is the charge density and B(r) is the magnetic field defined on sites of 
the dual lattice, i.e., a lattice obtained from the original (direct) one by translating 
its set of vertices to the centers of each plaquette of the direct lattice: 

B = eijAiAj (2.31) 
AiAj = Aj(x + et) - A(x) (2.32) 

The Gauss law constraint (2.30) implies that for each fermion on the site x of the 
direct lattice, there is also a quantum flux (or vortex) in the corresponding site r 
of the dual lattice. It can be implemented at the field theory level by coupling the 
fermions to an Abelian statistical Chern-Simons (CS) Gauge field. In order to show 
how to do it, let us consider the Lagrangian: 

L 2 d = ^ ip\x)iD ip(x) -t \^{x)e iAi i){x + e j ) + he ] 

x x,j=l,2 

+ \ J2e^, x A^(x)F^(x) , (2.33) 

x 

where 

D = d - iAo (2.34) 
Fij = A { Aj - AjAi (2.35) 
F 0l = d Ai - A t A . (2.36) 

The canonical quantization of the above Lagrangian in the Gauge Aq = imposes 
the constraint (2.30) at the level of the Hilbert space [17]. The classical solutions of 
the constraint can be written as follows: 

Aj =A 3 <P{x) = - e Y J [^ + e v r')-Q{x,r')]p{x'), (2.37) 

x' 
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where Q(x,r) is the lattice-angle function [18] which contains a branch cut from r 
to oo satisfying: j> r A6> = +1 for any closed curve T that encloses the point r of the 
dual lattice. Moreover, the operators: 

a(x) = e i * (x V(a:) (2.38) 
a\x) = i)\x)e- i ^ x) (2.39) 

satisfy bosonic commutations relations for 1/29 = mir and due to the Pauli principle, 
the operator (2.38) are hard-core bosons representing charge density waves (CDW) 
states of the underlying electrons. So that the mapping defined by the equations 
(2.38) (2.39), relate the fermionic degree of freedom to bosonic ones, and therefore 
the spinless fermion system can be mapped onto the following Hamiltonian: 

Hit = ~t [a\x)a(x + e„) + h.c] + U^p(x)p(x + a) . (2.40) 

x, fi=l,2 x 

The identifications: 

S+(x) = a\x) (2.41) 

S~(x) = a(x) (2.42) 

S z {x) =p{x) - 1/2 (2.43) 

where the operators S + (x) and S~(x) are the raising and lowering spin operators 
for a s = 1/2 spin particle, allows us to map the original fermion system onto the 
two-dimensional XXZ-spin model: 

Hfxz = E ["( + S " S j ) + AS * S j] ( 2 - 44 ) 

(ij) 

where i,j denote lattice sites and the sum is taken over nearest neighbors , and we 
have defined A = t/U 

Summarizing, fermionic two-dimensional systems can be bosonized by attaching 
fluxes to particles, which is achieved at quantum level by constraining the Hilbert 
space states after imposing the Gauss Law. The dynamics of these systems may 
be described at the Hamiltonian level in terms of these bosonic degrees of freedom, 
which physically represent charge density waves. Alternatively, the charge density 
waves may be replaced by another set of degrees of freedom, such as the spin waves 
in the so-called XX Z spin model. 



3 Integrability of the two-dimensional model 

As we have seen in the previous Section, the one- dimensional fermionic model with 
nearest neighbors interaction displays a Mott transition. This property can be estab- 
lished through the integrability of the (one-dimensional) XX Z spin chain, which is 
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equivalent to the linearity of the associated free boson system. Our strategy for dis- 
cussing the Mott transition in the corresponding two-dimensional spinless fermionic 
model on the square lattice (2.1) is to show that the property of integrability could 
also be extended to encompass this case. The discussion of the Mott transition in 
this system could then follow the line of reasoning of the previous Section. For 
reasons that will be clear latter, we will first consider the Six- vertex model [20]), 
and later show its relation to the lattice fermionic model. The Six-vertex model is a 
statistical model on a two-dimensional lattice, on which a classical electric current 
defined on each link can interact with other currents at the lattice sites. Each site 
of the lattice (which we will refer to as a vertex) may be in one of the six possible 
configurations shown in figure (1). The energy e v associated with a given vertex 
depends on the four current states at the edges only. If we further impose a parity 
(Z 2 ) invariance, we are left with three possible vertices: 

a = Wl = w 2 = e _/3ei (3.1) 
b = w 3 = w A = e' 1362 (3.2) 
c = w 5 = w e = e~ l3t:i (3.3) 

Following Refs [19] [20], one defines a vector space Vi for each vertical link, another 




wl w2 w3 

Figure 1: Boltzmann weights of the Six- vertex model 

one in a given horizontal row V a (the so-called auxiliary space) and a vertex operator 
R such that its matrix elements are interpreted as the Boltzmann weights of a given 
vertex according to: 

(jii+^OilRlfji^Pi) = W '(aiPi, fJ>i+ifJ>i) , (3.4) 

where |aj), denote the states on the vertical links, and the ones on 

the horizontal links. In this basis the R-Matrix of the Six-vertex model reads: 

" a 

b c 

c b 
_ a 

Following Baxter, we introduce the row-to-row transfer matrix which plays the role 
of a discrete-time evolution operator (the time variable is taken as flowing upwards 



(3.5) 
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from one row to the next). More precisely, the transfer matrix is an endomorphism: 

T : H N = V x <g> V 2 <g> ®V N ^V 1 ®V 2 ® ® V N , (3.6) 

defined by: 

T = Tr a [R aA rR a (jv-i) R a i] > (3-7) 

where the trace Tr a is taken on the auxiliary space. The integrability of the model 
is guaranteed by the existence of a set of mutual commuting row-to-row transfer 
matrices. In fact, it has been shown [20] that for two sets of Boltzmann weights, 
(a, b, c) and (a', b', c') the transfer matrices of the Six-vertex model satisfy: 

[T{a,b,c) , T(a',b',d) ] = & A 6v = A' 6v , (3.8) 

where A 6 „ = (a 2 + b 2 — c 2 ) / (2ab) (A' 6v is defined in an analogous way) is the so-called 
anisotropy parameter. This condition is equivalent to the existence of solutions of 
the Yang-Baxter equation [20]: 

Ri 2 (u)Ri 3 (« + v)R 23 (v) = R 2 3(v)'Ri3(u + v)R 12 (u) , (3.9) 

where u is the spectral (uniformization) parameter, and each operator Rj acts non- 
trivially on Vi<g>Vj. This equation can also be written in the form of the commutation 
relations for the quantum //-operators 

R 12 (u - v)L 13 (u)L 23 (v) = L 23 (u)L(v) 13 R 12 (u - v) . (3.10) 

As it is known [19], the transfer matrix of the Six- vertex model is related to the 
Hamiltonian of the XX Z model by: 

T(/i) = e -^-«( A ) , (3.11) 

which also shows that the Yang-Baxter equation implies the integrability of the 
XX Z spin model. 

Moreover, equation (3.9) may be considered as an equation among operators, 
whose solutions define integrable planar lattice models, i.e., (2 + 0) dimensional 
statistical systems or (1 + 1) dimensional quantum systems. Within this approach, 
the Six-vertex model may be considered as one specific solution of the Yang-Baxter 
equation when the space Vi is the representation space of spin 1/2 particles. Al- 
though a classification of the solutions of (3.9) is not known, some solutions have 
been found. These solutions are related to the quantum deformations of Lie alge- 
bras or, more precisely to the quantum deformations of the universal enveloping Lie 
Algebras also called quantum groups. In fact, the universal Yang-Baxter equation : 

R12R13R23 = R23R43R42 , (3-12) 
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which is the Yang-Baxter equation independent of the spectral parameter, has the 
so-called universal _R-matrices as solutions. For example, when V{ is the space rep- 
resentation of spin 1/2, Drinfeld [21] has given a solution: 



R 



H®H/2 



E 

n=0 



I -q~ )\^ q nH/2 {x+r ^ q -nH/2 {x _ r ^ 



(3.13) 



where q is a complex parameter, [n] q = (q n — q n )/(q — q x ) is a g-number and the 
generators X + , X_ and H satisfy the commutation relations: 



q 



H 



q 



-H 



q-q 



[H,X + }=2X, 
[H,X_] = 2X_ , 



(3.14) 

(3.15) 
(3.16) 



with co-products: 



A(X ± ) = X ± 
A(H) = H® 1 + 1 



q H/2 + q -H/2 



X. 



± 



H 



(3.17) 
(3.18) 



which define the quantum group U q (sl(2)) (note that the co-multiplication opera- 
tor A should not be confused with the anisotropy parameter, as is clear from the 
context). Furthermore, we define the operators [19] : 



E 1 



e"S- 



F 



Hi 



2S Z 
-2S Z 



(3.19) 
(3.20) 



x 



where , S' ± are the raising (lowering) operators of the spin-1/2 particle and 
is the afnnization parameter. These operators define an irreducible representation 

(e M , 1/2) of the afline algebra sl(2). In this context, the i?-Matrix act as an inter- 
twiner between the tensor product of two representations : 



R(e^,e^)A(g) = A'(g)R(e»\e»*) , 



(3.21) 



where g is any element of the Quantum Group and A' is the inverse co-product, i.e., 
the co-product composed with the operator permuting vector spaces. This i?-matrix 
has the form: 



R{e^\e^) 



qx — q 1 x 1 








x — x~ 
q-q- 






q-q~ 
x — x" 








qx — q _1 x~ 



(3.22) 



where x = e M , \x = H\ — \xi- This R- matrix coincides with that of the Six- vertex 
model for the parametrization a = sinh(w + ij) ,b = sinh(-u) c = i sin 7 q = exp (27) 
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(for details see [19]), so that the Six-vertex model possesses symmetry U q (sl(2)). 
Moreover, it have been shown in [22] that the Hamiltonian of the XXZ model in 

the thermodynamic limit commutes with the Affine Quantum Group U q (sl(2)) and 
that the space of states is identified with the tensor product of level 1 highest and 
level (—1) lowest representations of U q (sl(2)). Besides, the corresponding L-operator 
is a g-deformation of the fundamental L-operator of the XXX (Heisenberg) spin 
chain, given by [23]: 



T XXZ 

n,a 



CZ 1 _ c 

xq n — x q 1 



(q-q xq Sz - x l q s ™ 



(3.23) 



This Lax operator, together with the i?-matrix (3.22) satisfies the LLR 
condition (3.10)). It is possible to rewrite this condition introducing 



RLL 



L(x) = Q(x)L(x)Q-l(x) R = Q{x)Q(y)R(x/y)Q~\x)Q-\y) , (3.24) 



where 



Q(x) = 



x 



1/2 



X 



-1/2 



(3.25) 



which yields: 



R 



qx — q 1 x 1 










„-i 



x — x 
x(q-q~ l ) 






x~ 1 {q-q- 1 ) 
x — x^ 1 








qx — q~ 1 x~ 1 



(3.26) 



3.1 Three-dimensional structure of Quantum Groups and 
Vertex Models 

The (2 + l)-dimensional analogue of the Six-vertex model, is a 'quantum- vertex 
model' where the classical weights e _/3ea , e _/3efc e _/3ec , should be replaced by 'quantum 
vertex operators' defined on the lattice Hilbert space. In order to define this vertex 
model, we need first to consider a quantum lattice [24] defined as follows: for each 
lattice site (xi,yj), we define a Fock space Fij and a set of representation spaces of 
spin 1/2 particles V i: V}, V i+ i and Vj+i on each link joining two lattice sites. The 
states of the link in the lattice are arrows (as in the Six- vertex model), and the 
states in the Fock space are labeled by the number of particles in the site \riij) as 
shown in figure 2. Then, we assign a scattering amplitude (and a vertex operator) 
to each lattice site by: 

<s£S = W n v\^<y^\<*>PM) = L u» ( 3 - 27 ) 
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where Ly^y.^ is a 'three dimensional Lax operator' acting on the spaces Vi®Vj<8)F '. 
Furthermore, it is possible to define layer-to-layer transfer matrices T mn ({A}, {//}), 
where the pair (to, n) labels the rows and columns of a given layer and ({A}, {//}) 
are the spectral parameters, associated to the rows and columns, respectively, by: 

T mn = Try^yjJJ Y[L v .y. tF XXi, Hj)} = T mn {{X},{n}) , (3.28) 

where V x = ®™Vi and V y = <8>™V},{A} = {Ai, A 2 , ....A n }{//} = {//i, /i m }. Here, 
the layer-to-layer transfer matrix plays the role of a temporal evolution operator 
in an unitary time step, where the temporal axis coincides with the direction per- 
pendicular to the layer. For quantum systems in (2 + l)-dimensions (and also for 




Figure 2: Graphical representation of the sates \n = 0) and \n = 1) belonging to the Fock 
spaces at each lattice site, the arrows correspond to the states |oj) , |aj) ,\bj) ,\b'j) 

three-dimensional statistical systems), the integrability is guaranteed by the com- 
mutativity of the layer-to-layer transfer matrices. The integrability is also encoded 
in the the so-called Zamolodchikov Tetrahedron Equation (TE), which is the three- 
dimensional analogue of the Yang-Baxter equation, or equivalently in terms of the 
Lax operators, through the three-dimensional analogue of the LLR — RLL equation: 

RafecR-acieRfed/Rce/ = ^cef^bdf^adc^abc , (3.29) 
Ll2,aLl3 i f,L 2 3 lC R a f 1 c = Ra6cL23cLl3 iC Li3 i fe , (3.30) 
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where the operators Rijk define the mapping R abc '■ F a ® F b <E> F c — >■ F a <g> <8> i^j 
and the operators L act on V = 14 (8) V 2 <S> V 3 <g) F a <g> F b <g) F c . If the Fock space 
F a is considered as the representation space of some algebra A, then the operators 
L can be represented as operator-valued matrices acting on Vi <g> Vj such that their 
coefficients are given in terms of the generators v a of the algebra A and some complex 
parameters s a , so that the equation (3.30) takes the form of a local- Yang-Baxter- 
equation: 

£12 (v , s a )L 13 (v b , s 6 )L 23 (v c , s c ) = L 23 (v' c , s c )L 13 (v' b , s b )L la (v' a , s a ) . (3.31) 

For classical systems, the algebra A is chosen to be the Poisson Algebra P. In this 
case, the matrix L which solves the local Yang-Baxter equation (3.31) is given by 



where the indices in the second space enumerate the two-dimensional blocks , while 
those for the first space enumerate the elements inside the blocks and k — I — a* a . 
Moreover, 

{a*,aj}pB = 25ij {ki,dj}pB = Sijkidj \ki,a*}pB = —Sijkid* (3.33) 

where { , }pb denote the Poisson brackets. It was shown in [6] that equation (3.31) 
defines a canonical transformation (automorphism) of the triple tensor product of 
the Poisson algebra. This solution correspond to the classical three-wave problem, 
i.e., the linear propagation of three-dimensional waves along to three mutually per- 
pendicular axes. 

In quantum systems, we expect that the Algebra A will be the either bosonic or 
fermionic. This is actually the case for free Boson or Fermion systems. However, 
we are interested in interacting solutions of the TE. Recently, a new solution to the 
TE associated with the three-dimensional structure of the affine quantum group 
U q (sl(n)) has been found in [6] [7]. The new solution may be understood as the 
quantization of either the classical three-dimensional wave problem or the quantiza- 
tion of fluctuations of extended spatial objects. It amounts to taking L-operators as 
block matrices with two-dimensional blocks, in which matrix indices in the second 
space enumerate the blocks while those for the first space enumerate the elements 
inside the blocks. 



L lt 2(k a ,a a ,a* a ) 



10 

k a a* a 

-a a k a 

1 



(3.32) 



Lij(A v ) 




'V 









(3.34) 







-Q 



■v 
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with i,j — 1, 2, 3 , v — a, b, c and where now (k, a\ a) are quantum operators acting 
on the Hilbert space and the algebra A is the g-oscillator algebra defined by [6]: 



qa)a — q 1 aa) 



k 2 



a) a) 



q-q 

k = 



[h, a f ] 



,t 



[h, a] = —a 



(3.35) 
(3.36) 



The above solution allows one to define a 'quantum vertex model' by assigning 
'quantum vertex operators' fj at each vertex on the lattice, according the rules 
shown in the figure (3 ), where v 2 = —q~ 1 Xfi. We now map the square lattice onto a 
torus, and let C a and C& be the two basic homotopy cycles on that torus, such that 
each homology cycle corresponds to a path along one coordinate axis of the square 
lattice . Any path on the torus belongs to a given homotopy class P ~ nC a + mCf,, 
and it is possible to define [25] : 



T 

- 1 n. 



p jeP 



(3.37) 



where the sum over P means the sum over different paths, which is exactly the 
equation (3.28). The commutativity of the layer-to-layer transfer matrices follows 
from another (related) Tetrahedron Equation [26] [6]: 



M ii/ (H n/iJ l ')M jj/ (HoX/X')L ij (A v ,X,ij,)L i/r (A v ,X,ij l ) = 
Li'j'(A v , X,fi)L i:j (A v , X,ii)M jj ,(H X/X')M ii ,(H n/n') , 



(3.38) 



where the matrix elements of M(H ,() belong to an additional copy of the q- 
oscillator algebra denoted by H and : 



My (fib) 









X (-q() h » 













Vo( ho 



(3.39) 



with z/q = q ,_1 Ao/io- Equation (3.38) can be verified directly from the operator (3.39), 
and the commutativity of the layer-to-layer transfer matrices follows from its defini- 
tion and the use of equation (3.38). In the cubic lattice with boundary conditions( 
where the third dimension corresponds to the temporal axis) , the solutions of the 
TE given in (3.34) 3.39 have the following properties: 

• For q < 1, they yield the Fock space representation of the g-oscillator algebra: 



a|0) = \n) 



: |0) k\n) = q n+l/2 \n) 



(3.40) 



\fWvf) 

where (x; q 2 ) — (1 — x)(l — q 2 x) (1 — q 2 x). 

It is possible to construct the elements of the i?-matrix in the basis of the 
g-oscillator algebra (n^, n' 2 , n' 3 \R\ni, n,2, n%). The resulting i?-matrix is non- 
degenerate in F® . 
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V 2 



vK 



va* 



va 



Figure 3: Weights of the Quantum- Vertex-Model 

• The above solution implies that the standard Yang-Baxter equation is satisfied. 
One defines the operators 

Rbc = Tr Fa [R abl;Cl R ab2>C2 R a b n ,c n } (3-41) 

L v b = Try [Xo,i,&i-ko,2,&2 ^o,n,b„] (3.42) 

where F b = F bl ® F b2 ® ■ ■■■F bn , F c = F C1 ® F b2 ® ■■■■F Cn , which involves the 
product of n operators along the 'third direction, i.e., the direction labeled by 
'a'. Due to the fact that the i?-matrix is non-degenerate in F a ® F b ® F c , the 
TE implies that these operators satisfy 

RbcRbdRcd — RcdRbdRbc (3.43) 

LybLvcRbc — RbcLvcL v b ■ (3.44) 

This construction yields a projection from the TE to the YB ones. Using 
a similar projection, the equation 3.39 implies that there also exists another 
Yang-Baxter equation (for details see [6]): 

Rv^VjLvi^Lv^b = Lv itb L Vj , b Rvi,Vj ■ (3.45) 
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• The new solution is associated with the affine quantum group U q (sl(n)), where 
n is the range of the third dimension, i. e, n is the number of two-dimensional 
layers. 

For simplicity we now take A = 1 /i = 1 and the indices 1,2,3 denoting the quantum 
spaces a, b, c. Inserting the operator L (3.34) in the local Yang-Baxter equation 
(3.31), we obtain the explicit mapping i?i 2 3 given by: 

k' 2 a[ = k 3 ai — k\a 2 a\ 
a 2 = a 2 a 3 + k±k 3 a 2 (3.46) 
k' 2 a' 3 = k\a 3 — k 3 a\a 2 . 

Equations (3.46) are the Heisenberg equations of motion for the quantum operators 
(a,, a\, ki), where we have denoted with primes those operators that evolve forward 
in time t' — t + At which describe the time evolution of the quantum three-wave 
problem. 



k' 2 (a[)' = k 3 a\ — kia 2 a 3 
{a\)' = a\al + k\k 3 a\ 
k' 2 {a\)' = k\a\ — k 3 a\a\ 



3.2 Integrability and the two-dimensional Mott transition 

In this Section we will come back to the two-dimensional Fermion model defined by 
(2.1) to study its symmetries and integrability. We first recall (Section (1)) that it 
is equivalent to a two-dimensional XX Z spin model: 

h f = h xxz = - S?S x j + SfS y j - as?s; ] . 

ij 

On the other hand, the quantum vertex model defined in the previous Section has 
the remarkable property that can be projected from three to two dimensions. This 
means that the equations: 

LvbLvcRbc — RbcLvcLyb (3.47) 

LvifiLvjfiRvuVj = RviyjLvifiLv^b , (3.48) 

can be interpreted as arising from a two-dimensional system. Moreover, two remark- 
able properties of the quantum vertex model and the associated three-dimensional 
structure of quantum groups have been discussed in [6] : 

L Vb = ^ =l C{u k ,X,{fii}) (3.49) 
T m = Tr Vuk [C(w k , \m{fii}) C(u k ,\i,{fii})]. (3.50) 

The first equation shows that the operator Lyb decomposes into a direct sum of 
the fundamental ^-operators £ s '( n ) of the affine quantum group U q (sl(n)), where 
Uk is the highest weigh of the representation it Uk . The second one, shows that the 
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row-to-row transfer matrix of the quantum vertex model can be reconstructed from 
the fundamental /^-operators. 

The Quantum vertex model is stationary for L 12 (a, af,k) = L 12 (a f , a)', k'), i.e., 
for the case when the patterns for two consecutive time slices are identical so that 
the third dimension has range n — 2. In this case, the operator L v b is: 



u) 



and 



1 + «AiA 2 g fel+ft2 

C&fJl) 





Hi(q hl -u\ 1 \ 2 q ha - 1 ) 
-q~ l uX 2 ii 2 a\a 2 








^ 2 (q hl+h2 + q- 2 u\i\ 2 ) , 



(3.51) 



— q 1 X\^\a\a\ 
/i 2 (q h2 - uX!X 2 q hl 



(3.52) 



For hi — h 2 — 1/2 we have: 



q hl 



q 
1 



, q 



hi 



1 
q 



, aia\ 




1-q 2 



, a\a 2 



1-q 2 
1 



If Ai = A2 = 1 and H1 — H2 — 1, the last operator is the i?-matrix (3.26) of the XXZ 
spin chain( or Six-vertex i?-matrix) . Since is possible to project the quantum vertex 
model onto any lattice direction, any row or column should give rise to a XXZ spin 
chain, and one can expect that this system will be equivalent to a two dimensional 
XXZ model, (esta propiedad es algo asi como una generalizacion de la propiedad 
de invarianza modular) 



We are now ready to map among each other the original two-dimensional fermion 
model, the quantum vertex model and the two-dimensional XXZ spin model. The 
two-dimensional lattice fermion model describes charge density waves propagating 
on the lattice of the underlying electrons above ( and below) the half-filling state. 
Viewed at a fixed time, the wave vectors of these charge density waves cross them- 
selves at each lattice site, defining a 'Six-vertex model'. The half-filling condi- 
tion means that the Fock space of the g-oscillator algebra must be restricted to: 
{\n) = |1), I — 1}}. These states are in one-to-one correspondence with the states 
I t) I I) of the two-dimensional spin model. Therefore, the quantum ( Fock) space 
is also in correspondence with the representation space of a spin 1/2 particle. Intro- 
ducing this information in the equations (3.45) (3.51) (3.52), we see that for each line 
on the lattice we can define an XXZ (one dimensional) spin chain corresponding 
to the wave propagation of the CDW along this line. This fact allows us to identify 
the parameter A = U/t (representing the normalized Coulomb interaction) of the 
one-dimensional Fermion model with A = —(q + g _1 )/2, where q = exp (ry) is the 
deformation parameter of the quantum group. 
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Moreover, let us consider an interaction star , i.e., the interaction among a central 
spin, labeled by the index a an its nearest neighbors, labeled by the index j3. A 
direct calculation shows that: 

i d 

--^H 1 / a ( u ) Tr a,i3C{u) 1/2 (r a )C(u) 1/2 {r^] = -H a>p , (3.53) 

where a(u) = sinh(-u + ry), and we have taken the parameters u — v — 1. Tr a ^ is 
the trace over the space V ia <8> Vj a ® <8> and if Q/ g is the element (a, /3) of the 
two-dimensional XX Z Hamiltonian. Now taking : 

M*./*) = {l, t , , !/l ^-^ l>1 (3.54) 
I L(a, a 1 , k, r a ) r a ,rp are first neighbors 

we have that the partition function of the quantum vertex model (QVM) is: 

Z(QVR) ^Tr^]J[Tr v ^ Vp L(a a ,alk a ,r a )L(a ,alk,r )} , (3.55) 

a/3 

where, F nm = F® nn is the complete Fock space of the square lattice with n rows 
and m columns. Tracing over the vector space Vo (there is nothing special about 
the 'line zero', so that we can also trace over any row space Vi,V2....V n ) and using 
equation (3.42), we have: 

Z(QVM) = Tr Fnm l[[Tr Vfi (L vb (yp)L vb (yp +1 ) N )} (3.56) 



= Tr Fm H[Tr V e(R™ z (yp)R™ z (yp +1 )) N } (3.57) 

= Tr Fm [Y[e^ H ^] (3.58) 



= Tr F [e^ Ha ^} (3.59) 

= Z{Exxz) (3-60) 

where in the second line we have taken into account the fact that we are using the 

two-dimensional (spin 1/2) representation of U q (sl(2)). We arrive at the interesting 
identity: 

Z{QVM) = Z{H 2 x d xz ) , H 2 x d xz = H F . (3.61) 

This equation (3.61) shows that the two-dimensional XXZ spin lattice and the 
Lattice fermion system (2.1) have Quantum Group Symmetry U q (sl(2)). 

Summing up: the two-dimensional lattice Fermion model, viewed at two different 
times, can be considered as a two-layered three-dimensional system where the third 
direction coincides with the temporal axis (see figure 4). From an intuitive point of 
view, one can assign an arrow to the time propagation direction of the charge density 
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waves of the underlying fermions. At the stationary point, where the pattern of 
arrows does not change with time, these arrows define a vertex model with quantum 
group symmetry U q {sl{2)). Alternatively, one can say that due the projection-like 
property (3.43), the two-dimensional XX Z spin system may be decomposed in a 
consistent way into two one-dimensional chains, so that the entire system will have a 
Quantum Group symmetry ((U q (sl(2)))® N ). The property that (2 + l)-dimensional 
quantum systems in square lattice with periodic boundary conditions reduce to 
(1 + l)-dimensional quantum chains, implying that the {{U q {sl{2)))® N ) symmetry 
reduces to the U q (sl(2)) Quantum Group symmetry, was first noted in [6] [25]. 




Figure 4: Current patterns of the Fermion model at two different times. 
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4 Construction of the Effective Field Theory 



As we have shown in the previous Section, the two-dimensional fermionic system 

defined by (2.1) has a (2 + 1) Quantum-Group symmetry U q (sl(2)). Discrete time 
evolution is given by the transfer matrix T mn ({A„}, {/x m }). We have presented an 
explicit form for it in (3.28) which is a solution of the Zamolodchikov Tetrahedron 
Equation, rendering the model integrable in (2 + 1) dimensions. However, this solu- 
tion is not completely transparent from the physical point of view. To further study 
the system and this solution, we will follow the path of writing down an Effective 
Field Theory (EFT) formulation for it (for a review, see [3]), valid for long-distance 
and low-energy domains. EFTs provide useful frameworks for analyzing the behav- 
ior of many-body systems, specially near to a phase transition point, were universal 
properties are vastly dominating the physical properties of the systems under study. 
The general scheme for applying the method of EFTs starts by identifying the ef- 
fective degrees of freedom that dominate the low-energy regime of a given system 
and their characteristic symmetries. Note that these degrees of freedom are usually 
chosen by phenomenological reasons and could bear no resemblance with the micro- 
scopic degrees of freedom of any underlying model describing the system. One then 
proceeds to write down the most general local action in terms of second-quantized 
fields representing the selected degrees of freedom which are consistent with the 
noted symmetries. In the case at hand, the construction of the EFT is obtained by 
choosing fermionic fields as the degrees of freedom. By using the vertex model repre- 
sentation of the lattice fermion model (2.1), and its connection to the Chern-Simons 
theory (and analytical continuation in the Chern-Simons coupling constant given 
in [41]), we will now show that the corresponding Effective Field Theory is a dou- 
ble Chern-Simons theory with Gauge Group U(l) and Quantum Group symmetry 

U q (sl(2)) <g> U q (sl(2)), which is the symmetry of the exact solution (3.28) (3.34). 
4.1 The one-dimensional case 

For the sake of completeness and clarity of exposition, let us start by considering the 
EFT describing the one-dimensional Mott transition. As discussed in Section( 1), the 
one dimensional counterpart of the lattice fermion system (2.1) can be bosonized. Its 
effective degrees of freedom are the charge density waves of the underlying strongly 
correlated electrons, and for a general coupling (A) the system has a u(l) Kac- 
Moody symmetry [28]. The Mott transition point (A = 1) coincides with the 
L utter- Emery point, where the degrees of freedom are two chiral non- interacting 
fermions. At this point, the equation (2.20) can be written in coordinate space on 
a circle domain as: 

v f 27rR 

H= 2J dx [4( x )(- id x)Mx) + cc + ^(x^-id^x) + cc] , (4.1) 
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where ij) r {x) ( ipi(x)) are the right (left) relativistic chiral fermions (Weyl fermions) 
moving on a circle of radius R with speed v. Each chiral branch can be bosonized 
independently. Let us focus on one chiral component, say the right one. This can 
system has a quadratic action on a circle (of radius R = 1) [29]: 

7, pod />2n 

S = -1T dt (d t + vd x )<f>d x <f>, (4.2) 

4?r Joe J0 

where k is the coupling constant, (f>(x,t) is a real scalar field and x denotes the 
coordinate along the circle of length 2n and v is the 'speed of light'. By coupling this 
one-dimensional chiral system to an external electromagnetic field E, the equation 
of motion of the bosonic field changes to [30]: 

d t Q = ^E, (4.3) 

where E = d t A x — d x A t is the electric field pointing along the circle and A is the U(l) 
electromagnetic potential. Equation (4.3) displays the chiral anomaly of the Weyl 
Fermion. This description can be also obtained from the Abelian Chern-Simons 
theory with Lagrangian 

L cs = ^e^A,d u A p (4.4) 

defined on the disc D, whose boundary is the circle of radius R — 1. By making a 
Gauge transformation 5A p = ed^X, the action undergoes a variation: 

he f 

SS CS = — J \(d t A e - d e A t ) , (4.5) 

showing that Scs is n °t Gauge invariant on the boundary. Equivalently, there is a 
chiral current 

r = *k = > • (4 ' 6) 

that shows again the presence of the chiral anomaly of the Weyl Fermion for k — 1 
[30]. The chiral anomaly of one chiral component is canceled by the corresponding 
anomaly of the anti-chiral sector of the complete theory, which is anomaly free. 
This implies that the Mott transition in a one-dimensional system is described by a 
double Chern-Simons theory : 

S = A f (Fxe^'AfaAf - A f d z xe^ p A L ^d v A L p , (4.7) 

JDxR ' ' 47T J DxR 

which contains two chiral Gauge fields (A R , A L ) of opposite chiralities and where 
the value of the coupling constant is k — 1. 
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4.2 Two-dimensional Effective Field Theory 

Now let us consider the EFT for the Mott transition in the two-dimensional square 
lattice. As it has been remarked before, one must identify the correct degrees of 
freedom dominating the low-energy regime and their characteristic symmetries. The 
two-dimensional Jordan- Wigner transformation shows that the lattice fermionic sys- 
tem can be also considered as a bosonic system. We choose these bosonic degrees of 
freedom as the characteristic ones for the low-energy regime of the system. Further- 
more, these can be viewed in more physical terms as charge density waves, which we 
choose to construct the EFT. We look for an EFT with three-dimensional Quantum 
Group symmetry that is consistent with the properties established in Section ( 2.2). 
Namely, that it should posses a kind of vertex model interpretation, that it should be 
parity-invariant and that it must be projection-able onto one-dimensional theories, 
each giving rise to a one-dimensional Mott transition. We recall that in these one- 
dimensional theories the degrees of freedom split naturally in two non-interacting 
chiral branches when the interaction parameter becomes A = 1. We will now show 

that this EFT is a lattice double Chern-Simon theory with U q (sl(2)) ® U q (sl(2)) 
Quantum Group symmetry. 

4.2.1 Vertex Models and Chern-Simon Theory 

A direct connection between vertex models and Chern-Simons Gauge theory was first 
established in classical articles by Witten [41] [42] by using a non-Abelian Chern- 
Simons theory. Let us review this connection and consider a Gauge connection 
A = Aitf that belongs to an a Lie group G, where t\ are the generators of the 
group. The non-Abelian Chern-Simons theory is defined by the action: 



where M. is an orientated topological manifold, and trn denotes the trace on the 
representation R of the group G. The natural observables of this theory are Wilson 
loops: let us consider a link L as a disjoint union of the circles Ci and pick up a 
representation Ri for each circle. The expectation value of the Link can be calculated 

as: 



Here we can take Ai = HxR, where R represents the temporal axis and E is a 
Riemann surface. In [41] it has been shown that it is possible to define vertex 
models by replacing the classical currents living on the links of a lattice by Wilson 
lines in a Chern-Simons Gauge theory. Furthermore, it has been shown that the 
expectation values of this Wilson lines can be calculated from the data in the CS 




(4.8) 




(4.9) 
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theory. More specifically, taking the Gauge group 577(2) and projecting the three- 
dimensional knots onto the plane, the 'Boltzmann Weights' for these vertex models 
are given by: 

W up = q 1/2 5 Sl . S3fi - q^ 2 5 Sl+S2fi q (s2 - S3) (4.10) 
W und = <T 1/2 5 Sl - S3 ,o - q 1/2 5 Sl+S2fl q {s2 ~ S3) (4.11) 
W pc = e SuS2 q- s ^ 2 , (4.12) 

where W up (W un d) denotes the Boltzmann Weight for the vertex in which the line 
labeled by the representations of SU(2) denoted by s\ and S2, are above (below) the 
line labeled by s 3 and s 4 . Furthermore, W pc denotes the Boltzmann weights for the 
pair-creation, and q is the deformation parameter of the quantum group, which is re- 
lated to the coupling constant of the Chern-Simons theory by q = exp (in /k) . Note 
the a different factor in the definition of the deformation parameter of the Quantum 
Group with respect to [41]. Our definition match with the relations defining the 
Quantum Group (equations (3.14 ) (3.15)( 3.15) within the conventions adopted by 
Alvarez-Gaume et al. [47], while the definition of the deformation parameter in [41] 
agrees with the convention adopted in [45]. To make contact with the statistical 
Six-vertex model, one needs to take the Gauge group G = SU (2) and compute the 
'Bob amplitude'. It was show using Skein theory that any four coupling (including 
over-crossing , under-crossing and pair creation) is equivalent to the Bob-amplitude 

A = u.5 SuS3 5 S2tS4 + ve SltS2 e S3tS4 .q- {si+S3) (4.13) 

for some complex parameters u, v. Taking u = (qx' 1 — q~ 1 x) and v = (x — x' 1 ) the 
corresponding i?-matrix R Sl>Sj is given by: 

R(x, q) = (qx' 1 - g _1 :r)/ + (x - x~ 1 )U(q) (4.14) 
q 1 



U(q) = 



1 q- 1 



ifi<j (4.15) 



17(g) = 0, if i = 3 (4.16) 

which is a possible form of the i?-matrix for the Six- vertex model [45]. For a Gauge 

group G = SU(2), the Chern-Simons vertex models exhibit symmetry U q (sl(2)) in 
the same way that the classical Six- vertex model does. Therefore, the mathematical 
structure of the Quantum Groups encodes the topology of planar Wilson loops. A 
little of caution must be taken with the above defined Boltzmann Weights: first, note 
that since the Chern-Simons theory is well-defined for integer values of the coupling 
constant k, not all real values of the Boltzmann weights are well defined if we consider 
q = exp {jii/k). Second, the classical vertex model is defined in terms of classical 
degrees of freedom (the currents). However, the Chern-Simons action defines a 
non-dynamical theory (i.e., there are no dynamical currents in CS theory). Both 
problems are solved by the following argument: let us impose boundary conditions 
to the EFT, i.e., by compactifying the space domain onto a torus. Cutting down 
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the torus along any cycle induces a loose of the Gauge symmetry, so that the Gauge 
fields become dynamical degrees of freedom, and as it is well-known that the CS 
theory becomes then equivalent to a chiral Weiss-Zumino-Witten (WZW) model 
defined on a circle [42]) 

Scwzw — —k f d9dttr[g- l d e gg- l d Q g] + \ [ ^"Tr^d^gg-^gg-^g] , 

JCxR JT,xR 

, . (417) 

where g(z) is the chiral WZW field living in the group manifold of G. This is 
a conformal field theory (CFT) with central charge c — (k + \ SU(2) | ) /{k + c v ), 
where \SU(2)\ is the number of generators of the SU(2) Lie algebra, and c v is 
the dual coxeter number. The WZW model splits naturally into holomorphic and 
antiholomorphic pieces. The (holomorphic) energy-momentum tensor is given by 
the Sugawara form [43]: 

T{z) = * : J «(z)Ja{z) : , (4.18) 

^ v ' a 

where k is the Chern-Simons coupling constant and c v is the dual coxeter number, 
which for SU(2) is c v = 2. In other words, after quantization of the WZW model, 
there is a shift of the parameter k — > k + c v implied by the Sugawara construction. 
The currents J a {z) and the tensor T{z) may be expanded in modes in the usual way 
( for a review see [44]): 

J a (z)J>— V" (4.19) 

n 

T{z) = Y J Z- n - 2 Ln- (4.20) 

n 

The WZW model has conformal and affine su(2) k symmetries which can be written 
in terms of the Fourier modes of the Virasoro and current operators: 

[L n , L m ] = (n - m)L n+m + — 5 n+mfi (n 3 - n) (4.21) 

[Jl J b J = iV2e abc J c n+m + kn5 ab 5 n+mfi (4.22) 
[L n , J n a ] = -mJ a n+m . (4.23) 

The identification between the CS theory and WZW model was used in [41] to obtain 
an analytical continuation in k, which follows from the Knizhnik-Zamolodchikov 
equation: 



& ~ IT— E -^Mzi, zi)~~9{zN, z N ) = . (4.24) 
K + g Zi — Zj 



This can be used to compute the braiding matrices of this CFT. It has been shown 
that these braiding matrices correspond to the Boltzmann Weights of the associated 
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'interacting round a face '(IRF) model [47], which is equivalent to a vertex model 
(this is the so-called face- vertex equivalence). The Boltzmann Weights can now be 
defined by the equations (4.12), with the redefinition q = exp (m/(k + 2)). As it is 

known, the WZW model possesses a quantum group symmetry U q {sl{2))®U q {sl{2)) 
[47], with q = exp (m/(k + 2)) 

4.2.2 Abelian Chern-Simons Theory and Quantum Groups 

However, as it was pointed out by Witten [41], in the reduction from CS Gauge the- 
ory with Group G to any vertex model, one losses the local and Global G symmetry 
so that the vertex models retain only the maximal torus T symmetry of the Group 
G. For the case at hand where G = SU(2), the vertex model has a U(l) symme- 
try which naturally corresponds to the Gauge symmetry of the charged fermions 
(charge density waves) propagating on the lattice. The extended Kac-Moody sym- 
metry su(2) 1 present in the WZW model is obtained by taking the equivalent CFT 
of (two) chiral bosons at self-dual point i.e., where the currents: 



satisfy the algebra su(2)k ( equation (4.22)) at level k — 1. As we have seen in the 
previous Section, charge density waves are also naturally accounted by a Abelian 
Chern-Simons theory on the circle. 

Imposing periodic boundary conditions on the square lattice amounts to consider 
the Abelian CS theory on a (spatial) torus, such that the square lattice on the plane 
defines the homology cycles on the torus. In this domain, new degrees of freedom 
associated to the global Gauge transformations arise. To be more precise, let us 
consider the Lagrangian of the Abelian CS theory: 



where T denotes a torus with modular parameter r and homology cycle basis C a Cp. 
As it is known ([33] [32]), the Gauge field A can be parametrized in this domain by 
using the Hodge decomposition, which incorporates the windings around the non- 
contractible loops on the torus. In holomorphic coordinates z = x + iyz = x — iy 
it is given by 



(space independent ) function on time. The it can be shown that the Lagrangian, 
in the Gauge, A = becomes[32]: 



J z = id<j> , 



(4.25) 
(4.26) 




(4.27) 




a 




a 



(4.29) 
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The second term in the Lagrangian corresponds to Abelian CS theory on the plane 
with coupling constant k. The first one shows that the degree of freedom labeled 
by the function 'a' behaves as the coordinate of a quantum mechanical particle 
moving into an effective magnetic field B e ff = nkjlmgij) restricted to the lowest 
Landau Level. The quantization of this theory is well known and can be done in 
the Schodinger picture [31]. The wave functional is: 

^]=^(o) , (4-30) 

where ip(x) = e Xa - X e ~ /M, and X = y d± = (d 1 ±id 2 ). Small Gauge 

transformations do not affect the wave functional. However, due to the existence 
of non-contractible loops on the torus, the global Gauge transformations defined by 
[34] : 

a — > a + rii + rn 2 x ~ > X (4-31) 
A- z ->• -id-zRRT 1 (4.32) 

_. TliTY f Z 7127T f Z 

R(z,z) = exp{ / (uj — uj) / too — (too)) (4.33) 

Imr J Z0 Imr J ZQ 

affect both the zero and a modes. These large Gauge transformations are precisely 
the magnetic translations across a parallelogram unit cell: 

T R ij){a) = e*% laARl ^{a + R) , (4.34) 

where: 

T R = e ^ +ieA)R T a T b = e iB/2|aAb| T a+6 . (4.35) 

The exponential factor involves the flux through the parallelogram defined by a and 
b. Now, following [35] [36] we take the combination of the magnetic translations : 

E = 1 —-[T(a,a)-T(-a,a)] (4.36) 

Q-Q 

K = T(a,0) (4.38) 

where the translations are made from the elemental square plaquette of side a. These 
operators satisfy the relations: 

K - R- 1 

[E,F\ = -— r (4.39) 

kE = q 2 EK KF = q~ 2 FK , (4.40) 

that define the Quantum Group U q (sl(2). Here q = exp (iB/2a 2 ) = exp(i$/2) is the 
deformation parameter of the quantum group, and $ is the flux per plaquette. Note 
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that for $ = 2n we have q — — 1 and then 

T M) V^) = -ij>(x + a) (4.41) 
T (2a ,o)^(^) = i^(x + 2a) . (4.42) 

This is reminiscent to the staggered flux phase [37], where the fluxes are antifer- 
romagnetically ordered. Alternatively, the quantum group symmetry can be also 
analyzed directly at the level of the Gauge transformations acting on the wave 
functionals and the resulting quantum group is also U q {sl2) where the deformation 
parameter is identified directly in terms of the coupling constant k of the CS theory 
as q = e m l k [38]. Summarizing, the Gauge invariance of the Chern-Simons theory on 
the torus implies the existence of a Quantum group symmetry hidden in the theory. 
Now following [19] we can define the generators: 

E = e u E E 1 = e u E (4.43) 
F = e~ u E F 1 = e~ u E (AAA) 
K = K~ l K x = K , (4.45) 

where x = e u is an affinization parameter. These operators define a representation 
(e u , 1/2) of the affine Quantum Group U q (sl2). So that the double Abelian Chern- 
Simons Theory on the torus possesses a Quantum Group symmetry U q (sl2)<S>U q (sl2) , 
which is identified with the quantum group symmetry of the Six-vertex model or, 
equivalently, with the symmetry Group of the two-dimensional lattice fermion model. 
To write down the corresponding EFT on the square lattice with periodic boundary 
conditions all we need to do is to take the modular parameter r = % and to restrict 
the motion of the effective degrees of freedom (charge density waves) to the links 
on the lattice by replacing the Chern-Simons term with the Lattice Chern-Simons 
term ( which posses lattice differential operators). Therefore, the EFT is defined by 
the action: 

Sdcs = ^J d * x a l K w<% ~^J d3x a i K w a » > ( 4 - 46 ) 

with K^ v = S mu e^ u d a) S^f(x) = f(x + ae M ), d lx f(x) = (f(x + ae M ) - f(x))/a, 
(where a is the lattice spacing), which can also be written as a mixed Chern-Simons 
theory 



4.2.3 Identification of the Mott point (q = —1) 

We are now ready to identify the Mott point in the two-dimensional case. As we 
have seen in Section 3.2 the two dimensional Fermion model (2.1) which is equivalent 
to the XX Z spin lattice, can be split into one dimensional systems for each row or 
column of the square lattice , with Hamiltonian 

H l x d xz = + SfS^ - + H b , (4.47) 
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where Hi, = a{S{ — S Z L ) and a = (q — q 1 )/2. 

On the one hand the critical point of this system, which has been analyzed in 
Section (1), allows us to identify the Mott transition point in the square lattice with 
the one-dimensional transition at A = 1 . Since we have defined A = — (q+q~ 1 )/2, it 
implies q=-l. Besides the critical point in the one-dimensional system was identified 
with the WZW model or a Double Chern Simons theory on the circle. On the 
other hand we have identified the critical point of the two dimensional (double 
periodic) lattice fermion system (2.1)with the Double Chern Simon theory on the 

(space) torus with Quantum Group symmetry U q (sl(2)) <g>U q (sl(2)) . As we discussed 
in the Section (4.2.1), cutting down the torus the Gauge field become dynamical 
currents represented by a WZW model, which naturally splits its degrees of freedom 
in chiral components) with quantum group symmetry U q (sl(2)) ®U q (sl(2)) [47] with 
q = exp (m/{k + 2)). However, since the identification of the deformation parameter 
q is done at the classical (Hamiltonian) level, there is no Sugawara shift in the WZW 
coupling constant k due vacuum renormalization, an we should use q = exp (in/k). 
Therefore, since we can cut down the torus on any cycle, the affine Quantum group 
symmetry of the double Chern-Simons theory on the domain T x R implies the 
quantum Group symmetry of each XX Z spin system cut down in any circle. We 
therefore conclude that q = — 1 is the correct value of the deformation parameter 
characterizing the Mott point. 

The above statement is the EFT formulation of the statement discussed in [6] and 
reviewed in Section 3: the Quantum TE implies the standard Yang-Baxter equation, 
so that the integrability in the (2 + l)-dimensional quantum system implies the 
integrability of the (1 + l)-dimensional systems obtained from it by the projection- 
like character of the solution (3.34) . Sumarizing, the EFT of the lattice Fermion 
system (2.1) at the Mott transition point is represented by a double Chern-Simons 
theory at coupling constant k — 1, in agreement with the analysis based on the 
magnetic algebra presented in the Section 4.2.1. 

4.3 The Order Parameter 

We wold now like to discuss the emergence of the order parameter characterizing 
the transition. We can do this either from the point of view of the double-Chern- 
Simons theory or from the point of view of the equivalent CFT. First, we recall that, 
as we stated in the Section (3.2) , each row and column of the lattice system (2.1) 
can be considered as an one-dimensional system, which exhibit charge density wave 
order as we have stated in Section 1. At A = 1 + e the one dimensional system in 
the row is in gap-full state which correspond to the antiferromagnetic phase of the 
equivalent XX Z spin chain. For A = 1, the one-dimensional fermion system is in 

a gapless phase and is described by a CFT with su(2) 1 ® su(2) 1 symmetry, which 
is encoded in the WZW action or in a chiral- ant ichiral bosonic system compactified 
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on a circle at the self-dual radius. Let us now consider the transition in the square 
lattice: since we can cut along any cycle of the torus (defined in Subsection (4.2.2), 
we can obtain the above bosonic CFT at the self-dual radius in any row or column 
in the lattice. Invariance under surgery of the states in the CS theory, implies by 
consistency that the fermion system on the lattice should be described by a CFT 
with c—1. This is interpreted as the theory of the free boson on the (space) torus, 
with partition function: 

Z = _L_ q l/2(e/R+mR/2)-l/2{e/R-mR/2) ^ /^gx 
v ' e,m 

where rj{r) is the Dedekind function, r is the modular parameter of the torus and e 
and m are the electric and magnetic charges. Under i?-duality {i.e., the interchange 
R -H- 2/R) the partition function is invariant if one exchanges the electric and 
magnetic charges (e -h- m). 2 . This is also the partition function of an equivalent 
two-dimensional Coulomb gas described by the action: 

S ™ = \ + m ^)G{Rj ~ Rk)(^= + m ky fg)\ , (4.49) 

where G is the two-dimensional lattice Green function. Note that the above action 
describes the charge-charge , vortex-vortex and charge-vortex interactions. The 
Mott transition on the tours (square lattice with periodic boundary conditions) is 
given by the above partition function at the self-dual radius R = \/2, such that the 
system exhibits electric-magnetic duality. 

We would now like to discuss the behavior of the EFT away from the Mott 
critical point. We consider first the one- dimensional theory describing one of the 
orientations on the lattice. For A > 1 , the massive antiferromagnetic phase is 
represented in the continuous limit by a Sine Gordon Theory: 

S SG = J d 2 x[<^0<9 M + 2a o cos/30] . (4.50) 

It is well-known that the partition function may be written as [48]: 

Z = l } r §52 (^)]2 J nd 2 Zl exp(—J2m Ln \ z i ~ z i\ 2 + £ 2 ) ■ ( 4 - 51 ) 

Here z = x + iy denotes the position of the charges ^ in complex coordinates, that 
take values ±1. The renormalized coupling constant a = ao(e 2 )s^ plays the role 



2 By including all possible boundary conditions (periodic and anti-periodic) one obtains the 
S1/Z2 orbifold partition function that is also related to the Six- vertex model at the critical line 
[46]. We are not focusing in this case because we are working at half-filling states, which fixes the 
boundary conditions to be periodic 
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of a fugacity ( e is a short distance regulator). This describes the antiferroelectric 
phase of the Coulomb gas model. 

Moreover, the staggered flux phase revealed by the quantum group analysis of the 
previous Section can be described by a CS theory defined on a torus with punctures, 
by defining pseudo-spins on the dual lattice representing these fluxes. If we define 
the dual system as the two-dimensional XX Z model made from these pseudo-spins 
with coupling constant A' = A -1 , then the self-dual point is defined by A = 1. 
Here the self-duality in defined by invariance under exchange among the spins with 
coupling constant A in the direct lattice, and spins in the dual lattice with coupling 
constant A'. Just below the Mott transition ,ie., for A = 1 — e, the dual system is 
in the frozen (Mott) state with the pseudo-spins in an Neel state, so that the fluxes 
are also antiferromagnetically ordered. In this a case the EFT is given by a : 

v 

- 6(x d + l,y)-S(x d ,y d + l) + 5(x d + l,y d + l)} , (4.52) 

where a\ is an Abelian CS field and means that the sum is taken over all 
fundamental domains. Each domain has period 2a and contains four vortices in 
antiferromagnetic array. The emergence of only one CS term reflects the breakdown 
of the chiral symmetry (as can be seen comparing the equation (4.52) with the 
equation (4.46) which contain two chiral fields ), and the classical low-lying states 
reproduce an staggered -flux phase current pattern . This fact coincides with the 
quantum group analysis at q — — 1 presented in the Section (4.2). At the quantum 
level, Gauss law selects the physical states from the lattice CS Gauge theory on 
the torus with punctures. Therefore, the quantum order of the ground state of this 
theory is characterized as a staggered flux phase. 



The EFT that we have presented above allows one to study some properties of 
the system under doping. By analogy with the CS theory of the quantum Hall 
effect, we could expect a ground state stable against small doping. In that case, 
for the simplest inverse filling fractions k = m (m odd integer), the ground state is 
described as a droplet of incompressible quantum liquid [49] (however, other phases 
with more exotic quantum orders, like Nematic phases are also possible in other 
regimes ([52]) ([53]) ([54]) and is stable under small perturbations away from the 
center of a given plateau in the conductivity. In the Mott system, we have already 
assumed that the dynamically generated vortices act as external statistical fields for 
the new electrons injected in the system by doping (this can be considered as an 
extension of the /^-duality). At the self-dual point, statistical magnetic fields can be 
interchanged with statistical electric fields (on a torus). After imposing the lattice 
symmetries, the low-lying effective Hamiltonian for the injected electrons (in first 
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quantization) is: 



(4.53) 



where A, can take the values ±A. Therefore, the electric potential changes sign in 
x = ±y, producing domain walls between regions with different electron densities. 
Similar results can be obtained using the W4 symmetry, which is related to the 
relevant perturbations of the Ashkin- Teller and Six-vertex models away from the 
critical point [50] [51]. 

One consequence of having discussed the EFT is that, a-posteriori, the behavior 
of the electrons can be more easily understood. It can be shown that the interac- 
tion term in the Hamiltonian (2.1) in the the continuum limit contains a chemical 
potential term of the form —fj, p, with fj, — A , which ensures the half-filling con- 
dition. Therefore, changing the chemical potential by doping in 8/j, modifies the 
Hamiltonian (in the spin representation (2.6) by: 



For A = l,the dynamics of the electron system is given by the double CS theory 
(4.46), whose Hamiltonian can be defined as the temporal component of the stress- 
energy tensor H cs = T 00 , where T^ u = SScs/^Q^u and g^ u is the metric tensor. 
However, the CS action is topological and, therefore, independent of the metric im- 
plying H cs = for each chiral component, which leads to H(A = 1) = 0. This means 
that doping the system away from the critical point, the dynamics is controlled by 
an effective Ising Hamiltonian. 

5 Conclusions 

In this paper we have studied the Mott transition in an interacting electron system 
with a hopping term and nearest neighbors density-density coupling defined on a 
square lattice (2.1) . We have first reviewed the one-dimensional case with periodic 
boundary conditions ( i.e. when the system lives on a circle) in Section 1, and we 
have written down a Conformal Field Theory description leading to the identification 
of the Mott transition point (i.e., when the coupling constant is A = 1) as the 
Lutter-Emery point in the bosonic formulation, and the degrees of freedom with 
charge density waves. We have also pointed out that the low-energy dynamics 
at the transition is described by a Wess-Zumino-Witten model, as it was already 
implicit in the literature. 

We have also discussed the two-dimensional Mott transition, starting with the 
study of the integrability of the two-dimensional fermion system (2.1). To do that, 




(4.54) 
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we have used a two-dimensional Jordan- Wigner transformation and a new solution of 
the Zamolodchikov Tetrahedron Equation, which allowed us to identify the Fermion 
system (2.1) as a 'Quantum vertex model', and shown that the fermion system is 
is characterized by an affine Quantum Group symmetry. As a consequence of the 
projection-like property of the new solution, we concluded that the two-dimensional 
lattice system factorizes into two one-dimensional systems, one for any row and one 
for any column of the two-dimensional square lattice. This fact allowed for the 
identification of the two-dimensional Mott transition with the one-dimensional one, 
which occurs when the coupling parameter is A = —(q + q~ 1 )/2 = 1. The iden- 
tification of the symmetry and of the effective degrees of freedom (charge density 
waves) led to the construction of the Effective Field Theory at the critical point, 
whcih is a double (Abelian) Chern-Simons theory with Quantum Group symme- 
try U q (sl(2)) Cg> U q (sl(2)) and deformation parameter given by q = exp (— in/k). 
Furthermore, this effective theory may be considered as the broken phase of a non- 
Abelian Chern-Simons theory associated to the vertex models (which has been al- 
ready pointed out by Witten). 

Finally, the behavior of the system near the Mott point has also been investi- 
gated using the ideas of EFT. We have found that the transition is of the Kosterlitz- 
Thouless class, characterized by an array of Chern-Simons vortices in a anti-ferromagnetic 
order. This description corresponds to a d- density- wave order parameter for the 
matter currents. Upon doping with electrons, the magnetic-electric duality of the 
KT transition implies the appearance of domain walls between region of different 
densities. 
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